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We propose a new concept of modulated bipower variation for diffusion models with microstruc- 
ture noise. We show that this method provides simple estimates for such important quantities 
as integrated volatility or integrated quarticity. Under mild conditions the consistency of mod- 
ulated bipower variation is proven. Under further assumptions we prove stable convergence of 
our estimates with the optimal rate n^^''*. Moreover, we construct estimates which are robust 
to finite activity jumps. 
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1. Introduction 

Continuous time stochastic models represent a widely accepted class of processes in 
mathematical finance. Ito diffusions, which arc characterised by the equation 

Xt^Xo+ f a.ds+ /" asdWs, (1.1) 
Jo Jo 

are commonly used for modelling the dynamics of interest rates or stock prices. Here W 
denotes a Brownian motion, a is a locally bounded predictable drift function and a is a 
cadlag volatility process. A key issue in econometrics is the estimation (and forecasting) 
of the quadratic variation of X 

IV = f cr^ds, 
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which is known as integrated volatility or integrated variance in the econometric litera- 
ture. In recent years the availability of high frequency data on financial markets has moti- 
vated a huge number of publications devoted to measurement of the integrated volatility. 
A typical way to estimate the integrated volatility is to use the realised volatility (RV), 
which has been proposed by Andersen, BoUcrslcv, Dicbold and Labys [3] and Barndorff- 
Nielsen and Shephard [7]. RV is the sum of squared increments over non-overlapping 
intervals within a sampling period. The consistency result justifying this estimator is a 
simple consequence of the definition of the quadratic variation (see, e.g., Protter [21]). 
Theoretical and empirical properties of the realised volatility have been studied in nu- 
merous articles (see Jacod [17]; Jacod and Protter [19]; Andersen, BoUerslev, Diebold 
and Labys [3]; Barndorff-Nielsen and Shephard [7] among many others). 

More recently, the concept of realised bipowcr variation has built a nonparamctric 
framework for backing out several variational measures of volatility (see, e.g., Barndorff- 
Nielsen and Shephard [8] or Barndorff-Nielsen, Graversen, Jacod, Podolskij and Shephard 
[5]), which has led to a new development in econometrics. Realised bipowcr variation, 
which is defined by 

n-l 

BV{X,r,l),, = ^('^+0/2-1 ^ \A^xnA^^,X\\ (1.2) 

i=l 

with A"X = Xi/„ — X(j_i-)/„ and r, ^ > 0, provides a whole class of estimators for different 
(integrated) powers of volatility. Another important feature of realised bipower variation 
is its robustness to finite activity jumps when estimating the integrated volatility (in the 
case r V / < 2). This property has been used to construct tests for jumps (sec Barndorff- 
Nielsen and Shephard [9] or Christensen and Podolskij [11]). 

However, in finance it is widely accepted that the true price process is contaminated 
by microstructure effects, such as price discreteness or bid-ask spreads, among others. 
This invalidates the asymptotic properties of RV, and in the presence of microstructure 
noise RV is both biased and inconsistent (see Bandi and Russell [4] or Hansen and 
Lunde [15] among others). Nowadays there exist two concurrent methods of estimating 
the integrated volatility in the presence of i.i.d. noise. Zhang [24] has proposed to use a 
multiscale estimator as a generalisation of the concept of two scale estimators, which was 
introduced by Zhang, Mykland and Ait-Sahalia [25] based on a subsampling procedure. 
Another method is a realised kernel estimator, which has been proposed by Barndorff- 
Nielsen, Hansen, Lunde and Shephard [6]. Both methods provide consistent estimates of 
the integrated volatility in the presence of i.i.d. noise and achieve the optimal rate 
(whereas the two scale approach achieves the rate n~^/^). However, these procedures can 
not be generalised in an obvious way in order to obtain estimators of other (integrated) 
powers of volatility, such as the integrated quarticity, which is defined by 

IQ^ f atds. 
Jo 

This quantity is of particular interest because, properly scaled, it occurs as the conditional 
variance in central limit theorems for estimators of IV and has therefore to be estimated. 
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Moreover, both methods are not robust to jumps in the price process (here we would hke 
to mention the work by Fan and Wang [12], who obtain jump-robust estimates of IV by 
applying wavelet methods). 

In this paper we propose a new concept of modulated bipower variation (MBV) for 
diffusion models with (i.i.d.) microstructure noise. The novelty of this concept is twofold. 
First, this method provides a whole class of estimates for arbitrary integrated powers 
of volatility. Second, modulated multipower variation, which is a direct generalisation of 
MBV, turns out to be robust to finite activity jumps (when the powers are appropriately 
chosen). In particular, starting with MBV we construct estimators of IV and IQ which 
are robust to finite activity jumps. An easy implementation of MBV is another nice 
feature of our method. 

This paper is organised as follows: In Section 2 we state the basic notation and def- 
initions. In Section 3 we show the consistency of our estimators and prove a central 
limit theorem for their normalised versions with an optimal rate n"^/'*. In particular, we 
construct some new estimators of the integrated volatility and the integrated quarticity 
and present the corresponding asymptotic theory. Moreover, we demonstrate how the 
assumptions on the noise process can be relaxed. Section 4 illustrates the finite sample 
properties of our approach by means of a Monte Carlo study. Some conclusions and di- 
rections for future research arc highlighted in Section 5. Finally, we present the proofs in 
the Appendix. 



2. Basic notations and definitions 



We consider the process Y , observed at time points ti = i/n, i = 0, . . . , n. Y is defined on 
the filtered probability space (f2, J^, {J^t)t£[o,i]T P) ^^'^ exhibits a decomposition 

Yi/n^ X^fn + Ui, (2.1) 

where X is a diffusion process defined by (1.1), and {Ui)o<i<n is an i.i.d. noise process 
with 

Ea, = 0, EU^=uj^. (2.2) 

Further, we assume that X and U are independent. 

The core of our approach is the following class of statistics: 

M 

MBV{Y, r, 0„ = 1/2 ^ ) |^ |y Wj' ^ q, (2.3) 

771 — 1 

mn/M — K 
' i=[m-l)n/M 

with 

7^ m^/"^ 
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for some constants ci > and C2 > 1 (which will be chosen later). 

The intuition behind the quantity Y^j^^ can be explained as follows: Since X is a 
continuous process and (i7i)o<i<ra is an i.i.d. process with EUi = we deduce that 

^ mn/M-K 



n/M-K + l 

' i={m-l)n/M 
^ mn/M—K 



Y(^i+K)/n — ^{m-l)/M+K/n + Op{l). 



n/M -K + 1 

' i={m-l)n/M 

This means that the quantity Ym behaves like the increment X(m_i)/jv/_|_x/„ — 
X(^„i_iyM (although it has a bias that has to be corrected), and consequently it contains 
information about the volatility process a. 

The constants K and M control the stochastic order of the term Fm^'' . In particular, 
we have 




where U^f^ and are defined analogously to Fm^' in (2.4). By (2.5) the stochastic 

orders of the quantities in (2.6) are balanced, and we obtain 

Yi^^=Opin-'/^), (2.7) 

which explains the normalising factor in (2.3). 

More generally, we define the modulated multipower variation by setting 

M-k+l k 

MMnF,ri,...,r,)„ = n'-+/^"V2 ^ l[\Yi%]^,P , 

m=l j=l 

where fc is a fixed natural number, rj > for all j and r+ = ri + • • • + ffe . This type of con- 
struction has been intensively used in a pure Ito diffusion framework (see, e.g., Barndorff- 
Nielsen, Graversen, Jacod, Podolskij and Shephard [5] or Christensen and Podolskij [11] 
among others). Later on we will show that the modulated multipower variation for an 
appropriate choice of k and ri,...,rfc, turns out to be robust to finite activity jumps 
when estimating arbitrary powers of volatility. 

In the sequel we mainly focus on the asymptotic theory of the modulated bipower 
variation, but we also state the corresponding results for MMV{Y,ri, . . . ,rk)n for the 
sake of completeness. 
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3. Asymptotic theory 

In this section wc study the asymptotic behaviour of the class of estimators MBV{Y, r, Z)„, 
r,l>0. Before we state the main results of this section we introduce the following nota- 
tion: 

fir^E[\zn z^Af{0,l). (3.1) 



3.1. Consistency 

Theorem 1. Assume that £'|[/|2(''+')+'^ < oo for some e>0. If M and K satisfy (2.5) 
then the convergence in probability 

MBV{Y,r,l),,^MBV{Y,r,l) = ^ / {ly^al + ,^2i^^Y"^^^/^ du (3.2) 

ClC2 Jo 

holds. The constants vi and are given by 

ci(3c2-4+(2-C2)3v0) 2((C2-1)A1) 

3(C2-1)^ Ci(c2-1)^ 

Note that the limit MBV{Y,r,l) in (3.2) depends only on the second moment lu'^ of 
U, and no higher moments are involved. 

Next, we present the convergence in probability of the modulated multipower variation 
MMF(y,ri,...,rfc)„. 

Theorem 2. Assume that E\U\'^^+~^^ < oo for some e > 0. // M and K satisfy (2.5) 
then the convergence in probability 

MMV{Y, n, . . . , rk)n ^ MMV{Y, ri, . . . , r^) = f^IllZJ^ f\^^al + v^u^Y^I'' An 

ClC2 Jo 

(3.4) 

holds. 



3.1.1. Consistent estimates of integrated volatility and integrated quarticity 

Theorem 1 shows that MBV{Y,r,l)n is inconsistent when estimating arbitrary (inte- 
grated) powers of volatility. However, when r -I- / is an even number (this condition is 
satisfied for the most interesting cases), a slight modification of MBV{Y,r,l)n turns out 
to be consistent. Let us illustrate this procedure by providing consistent estimates for 
the integrated volatility and the integrated quarticity. 

As already mentioned in Zhang, Mykland and Ait-Sahalia [25] the statistic 

1 " 

'^' = ^El^V"-n-i)A/ (3-5) 
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is a consistent estimator of the quantity lo^ with the convergence rate n~^/'^. Conse- 
quently, we obtain the convergence in probability of the modulated realised volatility 

^1 Jo 

as a direct application of Theorem 1 and (3.5). 

Now we are in a position to construct a consistent estimator of the integrated quarticity. 
By (3.6) and Theorem 1 wc obtain consistency of the modulated realised quarticity, 
namely 

MRQjY).:^ (c.c./3)MW(r,4,0)„-2.,....^Mi^nn»-^l(^-)- T.^d.. 



(3.7) 

Note, however, that Theorem 1 provides a whole class of new estimators of the integrated 
volatility and the integrated quarticity. 



3.1.2. Robustness to finite activity jumps 

As already mentioned in the introduction, one of our main goals is finding consistent 
estimates of volatility functionals when both microstructure noise and jumps are present. 
For this purpose we consider the model 

Z = Y + J, (3.8) 

where y is a noisy diffusion process defined by (2.1) and J denotes a finite activity jump 
process, that is, J exhibits finitely many jumps on compact intervals. Typical examples 
of a finite activity jump process arc compound Poisson processes. 

The next result gives conditions on ri , . . . , under which the modulated multipower 
variation MMV{Z, ri, . . . , rfc)„ is robust to finite activity jumps. 



Proposition 3. // the assumptions of Theorem 2 are satisfied, max(ri, . . . ,rk) < 2 and 
Z is of the form (3.8) then we have 

MMV{Z, n, . . . , rk)n ^ MMV{Y, r,,..., r^), (3.9) 
where MMT^(y, ri, . . . , r^) is given by (3.4). 

Proposition 3 is shown by the same methods as the corresponding result in the noiseless 
model (i.e., U = 0). We refer to Barndorff-Niclscn, Shcphard and Winkcl [10] (Proposition 
1, page 799) for a detailed proof. 

Now we can construct consistent estimates for the integrated volatility and the inte- 
grated quarticity, which are robust to noise and finite activity jumps. Since o)^ is robust 
to jumps, the convergence in probability 

Mi^n^)..:-^'^^^^^"^^'^^^^^'^'^^""^^"'^ f'^l<iu (3.10) 

J^i Jo 
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holds as a direct consequence of Proposition 3. Similar to the previous subsection, a 
robust (tripower) estimate of the integrated quarticity is given by 

MTQ{Z)n 

(ciC2//x3/3)MMl/(Z,4/3,4/3,4/3)„ - 2i^iiy2U}^MBViZ)„ - i^U^^ 



(3.11) 



Remark 1. Recall that the realised volatihty RV converges in probability to the in- 
tegrated volatility plus the sum of squared jumps in the jump-diffusion model. It is 
interesting to see that the presence of jumps destroys the consistency of the estimator 
MRV {Z)m which can be interpreted as an analogue of RV . This is explained by the fact 
that jumps appear with different factors in the statistic MRV {Z)n, according to their 
positions in the intervals Y^^jf^, f^]- 

In contrast to our approach, the multiscale estimator of Zhang [24] and the realised 
kernel estimator of Barndorff-Nielsen, Hansen, Lunde and Shephard [6] converge in prob- 
ability to the quadratic variation of the jump-diffusion process X -\- J (in the presence of 
noise). 

Another important object of study is the impact of infinite activity jumps on the mod- 
ulated bipowcr (multipower) variation. Such studies can be found in Barndorff-Nielsen, 
Shephard and Winkel [10], Woerner [23] and in a recent paper of Ait-Sahalia and Jacod 
[1] for the noiseless framework. We are convinced that similar results hold also for the 
noisy model, although a more detailed analysis is required. 

3.1.3. Relaxing the assumptions on the noise process U 

So far we have assumed that U is an i.i.d. sequence and is independent of the diffusion X . 
Hansen and Lunde [15] have reported that both assumptions are somewhat unrealistic 
for ultra-high-frequency data. In the following we demonstrate how these conditions can 
be relaxed. 

First, note that the i.i.d. assumption is not essential to guarantee the stochastic order 
of Uln^ in (2.6). When we assume, for instance, that J7 is a g-dependent sequence, the 
result of Theorem 1 holds, although higher order autocorrelations of U appear in the 
limit. In this case we require a stationarity condition on U for the estimation of the 
autocorrelations and a bias correction of the limit in (3.2). 

Further, by using other constants M and K the influence of the noise process U can 
be made negligible, and independence between X and U is not required. (2.6) implies 
that, in particular, when we set 



K^cin^'^+\ M= — 
C2K 



(3.12) 
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for some < 7 < ^, the quantity Xm driven by the diffusion process dominates the 
term U^K More precisely, the convergence in probabihty 

^(i-27)(.+o/4-(i-27)/2y-|y^W|.|ym|;^Av^ / i^jr+^du (3.13) 

holds. The convergence in (3.13) has another useful side effect. It provides consistent 
estimates for arbitrary integrated powers of volatility. However, since the diffusion term 
Xjn^ dominates the noise term the above choice of K and M leads to a slower 

rate of convergence. 



3.2. Central limit theorems 



In this subsection we present the central limit theorems for a normalised version of 
MBV {Y,r,l)n- For this purpose we need a structural assumption on the process a. 
(V): The volatility function a satisfies the equation 

(Tt = ao+ f a; d.s + / a',dWs+ I v', dV^. (3.14) 
Jo Jo Jo 

Here a', a' and v' are adapted cadlag processes, with a' also being predictable and locally 
bounded, and F is a second Brownian motion independent of W. 

Condition (V) is a standard assumption that is required for the proof of the central limit 
theorem for the pure diffusion part X (see, for example, Barndorff-Nielsen, Graversen, 
Jacod, Podolskij and Shephard [5] or Christensen and Podolskij [11]). 

For technical reasons we require a further structural assumption on the noise process 
U. We assume that the filtered probability space {n,J^,{J^t)telo,i],P) supports another 
Brownian motion B ~ (5t)te[o,i] that is independent of the diffusion process X, such 
that the representation 

U^ = x/?itj(Bj/„ - B(j_i)/„) (3.15) 

holds. 



Remark 2. Condition (3.15) ensures that both processes X and U are measurable with 
respect to the same type of filtration. This assumption enables us to use standard central 
limit theorems for high frequency observations (see Jacod and Shiryaev [20]). The same 
assumption has already been used in Glotcr and Jacod [13, 14]. 

The normal distribution of the noise induced by (3.15) is not crucial for our asymptotic 
theory, and other hmctions of rescaled increments of B can be considered. Of course, this 
leads to a slight modification of the central limit theorems presented below. 

In the central limit theorems that will be demonstrated below we use the concept of 
stable convergence of random variables. Let us shortly recall the definition. A sequence 
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of random variables G„ converges stably in law with limit G (throughout this paper we 

write Gn — ^ G), defined on an appropriate extension {Q' ,J-' ^ P') of a probability space 
(Q,J^,P), if and only if for any JF-measurable and bounded random variable H and any 
bounded and continuous function g the convergence 

lim E[Hg{G„)] = E[Hg{G)] 

n — 'oo 

holds. This is obviously a slightly stronger mode of convergence than convergence in law 
(see Renyi [22] or Aldous and Eagleson [2] for more details on stable convergence). 
Now wc present a central limit theorem for the statistic MBV{Y,r,l)n. 

Theorem 4. Assume that U is of the form (3.15) and condition (V) is satisfied. If M 
and K satisfy (2.5), and 

1. r,le (l,oo)U{0} or 

2. r or I (0, 1], and ct, ^ for all s, 

then we have 

n^'^{MBV{Y, r, l)n - MBV{Y, r, I)) ^ L{r, I), 
where L{r, I) is given by 



V C1C2 Jo 

Here W denotes another Brownian motion defined on an extension of the filtered prob- 
ability space (ri, JF, (jFt)tg[o,i] , P), which is independent of the a-fi,eld T. 

Since o)^ — o;^ = Op(n~^/^) we obtain the central limit theorems for the estimates 
MRV{Y)n and MBV{Y)n defined by (3.6) and (3.10), respectively, as a direct conse- 
quence of Theorem 4. 

Corollary 1. Assume that U is of the form (3.15) and condition (V) is satisfied. If M 
and K satisfy (2.5) then we have 

n'/^(^MRV{Y)r.- I^ldu^ ^ l\,.,al + ^,cj')dWi, (3.17) 

where W' is another Brownian motion defined on an extension of the filtered probability 
space (fi, J^, (J^t)tg[o.i],P), which is independent of the a-fieldj-. 

Corollary 2. Assume that U is of the form (3.15) and condition (V) is satisfied. If M 
and K satisfy (2.5), and Us 7^ for all s, then we have 

n^/^(^MBV{Y)n- (Jldu^ 
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(3.18) 



4 9 



where W' is another Brownian motion defined on an extension of the filtered probability 
space (fi, JF, (J^t)tg[o.i],P), which is independent of the a-fieldj-. 

Now let us demonstrate how Corollaries 1 and 2 can be applied in order to obtain 
confidence intervals for the integrated volatihty. Note that the central hmit theorem in 
(3.17) is not feasible yet. Nevertheless, we can easily obtain a feasible version of Corollary 
1. Since the Brownian motion W is independent of the volatility process a, the limit 
defined by (3.17) has a mixed normal distribution with conditional variance 

vi Jo 

By an application of Theorem 1 the statistic 



Pl = ^MBV{YA,0)r 



is a consistent estimator of . 

Now we exploit the properties of stable convergence to obtain a standard central limit 
theorem 

^ ^—f^ — ^-^—^ --^JV{0,1). (3.19) 

From the latter confidence intervals for the integrated volatility can be derived. A feasible 
version of Corollary 2 can be obtained similarly. 

With the above formulae for and in hand we can choose the constants ci and 
C2 that minimise the conditional variance. In order to compare our asymptotic variance 
with the corresponding results of other methods we assume that the volatility process a 
is constant. In that case the conditional variance is minimised by 



y (c2 - 1)(4 - C2) cr 5 
and is equal to -2p= . (j3^ ~ 20.11a^uj. Note that the limits in Corollaries 1 and 2 arc the 

^ 3^18 

same up to a constant. Consequently, the asymptotic conditional variance of MBV (Y)n 
is minimised for the same choice of ci and C2 as above, and is approximately equal to 
26.14(T'^a;, when the volatility function is constant. 

As already mentioned in Ait-Sahalia, Mykland and Zhang [25] (see also Gloter and 
Jacod [13, 14]) the maximum likelihood estimator (when U is normally distributed) con- 
verges at the rate n~^^^ and has an asymptotic variance Sa'^w, which is a natural lower 
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bound. In contrast to our concept, the family of modified Tukey-Hanning kernel estima- 
tors as proposed by Barndorff-Niclsen, Hansen, Lunde and Shcphard [6] has an optimal 
asymptotic variance of about S.Oltr'^a;. This shows that our estimator is somewhat in- 
efficient in comparison to these kernel-based estimators. A natural direction of future 
research is to modify our procedure in order to achieve a higher efficiency. 

However, the concept of modulated bipower (multipower) variation has been estab- 
lished to provide estimates of arbitrary powers of volatility for the noisy diffusion model, 
which are additionally robust to finite activity jumps. These are properties which are not 
captured by the multiscale or realised kernel approach. 

For the sake of completeness we state a central limit theorem for the modulated mul- 
tipower variation MMV{Y,ri, . . . ,rfe)„. 

Theorem 5. Assume that U is of the form (3.15) and condition (V) is satisfied. If M 
and K satisfy (2.5), and 

1. ri,...,rfc e (l,oo)U{0} or 

2. r.i G (0, 1] for at least one i, and CTs 7^ for all s, 
then we have 

n^'\MMV{Y,ri,. . . , rfc)„ - MMV{Y, n, . . . , r^)) ^ L^n,. . . , r^), 
where . . . , r^) is given by 



L{rJ)^ 




(3.21) 



with 

k k fc— 1 j k k—j 

A{rl,...,rk)^Yl^l2n-i2k-l)Y[^4^+2^Yl^ir^ l^nJl f^r,+n+,- 

1=1 1=1 j=i 1=1 i=k-j+i 1=1 

Here W denotes another Brownian motion defined on an extension of the filtered prob- 
ability space {^l,!F ,{J-t)ti^[Q^i]T P) , which is independent of the a-field !F. 

Note that the constant A{ri,...,rk) also appears in the central limit theorem for 
multipower variation in a pure diffusion framework (sec Barndorff-Niclsen, Gravcrsen, 
Jacod, Podolskij and Shcphard [5]). 

4. Simulation study 

In this section, we inspect the finite sample properties of various proposed estimators for 
both integrated volatility and quarticity through Monte Carlo experiments. Moreover, 
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we compare our estimators' behaviour with the properties of the corresponding kernel- 
based estimators from Barndorff-Nielsen, Hansen, Lunde and Shephard [6]. To this end, 
we choose the same stochastic volatihty model as in their work, namely 

dXt^ lidt + atAWt, crt = exp(/3o + An), 

(4.1) 

drt art At + dBt, corr(dWt, dBf ) = p, 

with fi = 0.03, (3o = 0.3125, /3i = 0.125, a = -0.025 and p = -0.3. {U^)o<^<n is assumed to 
be i.i.d. normally distributed with variance w^. 



4.1. Simulation design 

We create 20 000 repetitions of the system in equation (4.1), for which we use an Eulcr 
approximation and different values of n. Whenever we have to estimate uj'^, we choose 
(D^ as defined in (3.5). 

Since we state propositions for a whole class of estimators, we do not focus on one 
special estimator. To be precise, we investigate the finite sample properties in three 
different situations. 

First, we study the performance of MRV{Y)n as an estimator for the integrated volatil- 
ity and compare it with the corresponding kernel-based statistic from Barndorff-Nielsen, 
Hansen, Lunde and Shephard [6], using the modified Tukey-Hanning kernel with p = 2. 
We denote this estimator by KB{Y)n- In Table 1 we present the Monte Carlo results for 
both mean and variance of the two statistics for various choices of n and = 0.01, 0.001, 
which is a reasonable choice, since IV is about 2 in model (4.1). Figure 1 gives histograms 
both of the standardised statistic in (3.19) and of the statistic 

nV4(log(Mi?F(y)„)-log(/>gdu)) r> 

PnlMRV{Y),, ^A/(0,1), (4.2) 

which is obtained via an application of the delta method. Both statistics converge weakly 
to a standard normal distribution. In this case, we choose two different values of n, namely 
n = 1024 and n = 16384. 

Second, we analyse the performance of the estimation of the integrated volatility in the 
presence of finite activity jumps. In this case we use the bipowcr estimator MBV {Z)m 
which is robust to jumps. We present the Monte Carlo results for this estimator in Table 2. 

At last, we analyse how well MRQ{Y)n works as an estimator for the integrated 
quarticity in contrast to the proposed bipower variation estimator in Barndorff-Nielsen, 
Hansen, Lunde and Shephard [6], which we call BP{Y)n- Note that BP{Y)n has a conver- 
gence rate of ri~^/^, which is obviously slower than the convergence rate of our estimator 
MRQ{Y)n- Table 3 shows the results in the quite simple setting 

AXt^ pdt + dWt (4.3) 



with jjL = 0.03 as above. 
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As mentioned before, the asymptotic (conditional) variance of the estimators MRV{Y)n 
and MBV {Y)n can be minimised for an appropriate choice of ci and C2, if the volatihty 
function is constant and the drift function is zero. However, even in model (4.1) and 
(4.3) the choice of ci and C2 as in (3.20) and with a replaced by IV may give an idea 
of a reasonable size for the two constants. Using a first estimate for IV we decided to 
choose ci = 0.25 for cj^ = 0.01 and ci = 0.125 for = 0.001, whereas C2 = 2. Since the 
computation of the optimal values of ci and C2 for the estimation of IQ involves the 
solution of polynomial equations with higher degrees than two, we have dispensed with 
this analysis and set ci = 1 and C2 = 1.6, both for = 0.01 and = 0.001. To produce 
the process J we allocate one jump in the interval [0,1]. The arrival time of this jump 
is considered to be uniformly distributed, whereas the jump size is Af(0, h?) distributed 
with /i = 0.1, 0.25. 

Note finally that it might be convenient not to plug in vi to compute MRV{Y)n and 
the estimator for the conditional variance 0f^^ but to use 

(C2-1)V" 

which is a better approximation to the second moment of n than vi. Similarly, 

one could use a finite sample analogue 1^2 ^2 as well. It has turned out that the 

performance of our estimators is indeed sensitive to the choice of vi (which is why we 
used J^i""*), whereas it is almost unaffected by the transition from 1^2 to 1^3"^ (which is 
why we used 1/2). 



4.2. Results 

Since our aim is mainly to give an idea of how well the different estimators work, we 
content ourselves with computing the estimated mean and variance of the bias-corrected 
statistics. Except for MRV{Y)n we therefore do not evaluate the accuracy of the stated 
central limit theorems. 

Table 1 shows that MRV{Y)n works quite well as an estimator of the integrated 
volatility in the noisy diffusion setting, since both bias and variance are rather small, 
at least for sample sizes larger than n = 1024. For large values of n and = 0.01 it 
provides even better finite sample properties than KB(Y)n, whereas the kernel-based 
estimator improves a lot when the variance of the noise terms becomes smaller. Never- 
theless, MRV{Y)n is a serious alternative to the kernel-based estimator, especially for 
large values of . 

Figure 1 indicates that the behaviour of each standardised statistic depends heavily 
on the actual size of w^, especially for a small sample size n. In fact, even for n = 1024 
we see a reasonable approximation of the standard normal density when the variance of 
the noise variables is large, whereas for ^ 0.001 the histo gram exhibits a significant 
shift to the right. One might suggest that these effects are caused by a large variance 
of the estimator of the integrated quarticity. For a large value of n all statistics work 
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Figure 1. Various histograms of the statistics defined in (3.f9) and (4.2). fn each line, the first 
data set was computed with n = f024, whereas for the second one we used n = f6 384. The first 
four histograms illustrate the case (3.f9) with — 0.01 (in the first line) and = 0.001 (in the 
second line) . The latter four graphics show the finite sample properties for the weak convergence 
in (4.2), in the same order as above. For comparison, the dashed line shows the graph of the 
standard normal density and the solid line gives a standard kernel density estimate. 
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Table 1. 


The Monte Carlo results for mean and variance 


of both MRV(Y)„ 


— J du and 


KB{Y)n ~ 


Jq '^u d« for various values of n and uj^ . The data 


are generated from the model (4.1) 




tj^ = 0.01 




2 

UJ = 


0.001 


71 


Mean 


Variance 


Mean 


Variance 


MRV{Y)^ 










256 


0.1363 


0.63 


0.5245 


1.782 


1024 


0.0433 


0.219 


0.1717 


0.269 


4096 


0.0113 


0.102 


0.0478 


0.055 


9216 


0.0045 


0.064 


0.0243 


0.031 


16 384 


0.0059 


0.05 


0.0129 


0.021 


25 600 


0.004 


0.039 


0.0094 


0.017 


KB{Y)„ 










256 


-0.022 


0.228 


-0.0289 


0.143 


1024 


0.0074 


0.091 


-0.0075 


0.042 


4096 


0.0195 


0.046 


0.0004 


0.015 


9216 


0.0203 


0.038 


0.001 


0.009 


16 384 


0.0201 


0.04 


0.001 


0.007 


25 600 


0.0178 


0.046 


0.0013 


0.005 



Table 2. Mean and variance of MBV{Z)n — du in the presence of jumps. We choose the 

sample frequency as before and analyse the finite sample properties for different values of lj^ 
and h, where h denotes the variance of the jump size 



n 


2 

LU — 


0.01, ^ = 0.25 


(.J =0.001, h 


= 0.25 


a; =0.001, h 


= 0.1 


Mean 


Variance 


Mean Variance 


Mean Variance 


256 


0.0582 


0.614 


-0.0839 


0.332 


-0.1224 


0.29 


1024 


0.0835 


0.295 


0.0274 


0.133 


-0.102 


0.112 


4096 


0.0707 


0.15 


0.0466 


0.063 


0.0184 


0.056 


9216 


0.0642 


0.102 


0.0461 


0.043 


0.0107 


0.038 


16 384 


0.0599 


0.076 


0.044 


0.032 


0.025 


0.028 


25 600 


0.0566 


0.059 


0.0415 


0.025 


0.0181 


0.023 



pretty well; however, it is remarkable that the transition to the log-transformed statistic 
provides an improvement, at least for a large choice of w^. 

From Table 2 we conclude that in the noisy jump-diffusion framework the proposed 
bipower estimator MBV{Z)n works quite well. Both bias and the variance of MBV{Z)n 
are rather small, even for moderate values of n. 

Finally, we see from Table 3 that MRQ{Y)n is on average closer to the true inte- 
grated quarticity than BP{Y)n, whereas the variance of BP{Y)n is smaller than that 
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of MRQ{Y)n, even though BP{Y)n has a slower rate of convergence. However, we are 
convinced that the efficiency of MRQ{Y)n can be improved massively by choosing the 
constants ci and C2 optimally. 

5. Conclusions and directions for future research 

This paper highlights the potential of the modulated bipower approach, and we are con- 
vinced that many unsolved problems in a noisy (jump-) diffusion framework can be tack- 
led by our methods. Let us mention some most important directions for future research. 
First, we intend to modify our approach by subsampling the statistic MBV{Y,r,l)n to 
obtain more efficient estimators of the integrated volatility and the integrated quarticity. 
A first step in this direction has been made in a recent paper by Jacod, Li, Mykland, 
Podolskij and Vetter [18] who proposed a subsampled version of MBV{Y,2,0)n to esti- 
mate IV in the presence of a more general noise process. Second, we plan to derive a 
multivariate version of the current approach. This can be used to estimate the quadratic 
covariation, which is a key concept in econometrics, in the presence of noise. Finally, 
an interesting and very important modification of this problem is the estimation of the 
quadratic covariation for non-synchronously observed data in the presence of noise (see 
Hayashi and Yoshida [16] for more details in a pure diffusion case). 

Table 3. The finite sample properties of MRQ{Y)n — a^du and BP{Y)„ — a^du in 
model (4.3). Both sample frequency and noise are the same as in Table 1 



Lj^^OM a;^ = 0.001 



n 


Mean 


Variance 


Mean 


Variance 


MRQ{Y)^ 










256 


0.0745 


1.348 


0.0686 


1.274 


1024 


0.0128 


0.587 


0.0121 


0.557 


4096 


0.0135 


0.306 


0.0013 


0.278 


9216 


0.0113 


0.203 


0.015 


0.184 


16 384 


0.0159 


0.152 


0.0155 


0.14 


25 600 


0.0088 


0.117 


0.0077 


0.108 


BP(y)n 










256 


-0.2517 


0.304 


-0.2803 


0.274 


1024 


-0.1811 


0.186 


-0.1434 


0.169 


4096 


-0.0312 


0.108 


-0.0745 


0.095 


9216 


-0.0089 


0.077 


-0.04 


0.065 


16 384 


0.0078 


0.059 


-0.0287 


0.048 


25 600 


0.0148 


0.047 


-0.0206 


0.039 
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In the following wc assume without loss of generality that a, cr, a', cr' and v' are bounded 
(for details see BarndorfF-Nielsen, Graversen, Jacod, Podolskij and Shephard [5]). More- 
over, the constants that appear in the proofs are all denoted by C (even if they depend 
on the actual choice of the powers I and r). 

Note first that the numbers vi and i>2 as given in (3.3) are the limits of the variances of 
the random variables n and 'n}^^U^\ respectively. Since a standard calcuation 

(n) 

shows that the true variances u!> are of the form 

i^i")=M, + 0(n-i/2), 

(n) 

one may conclude from the mean value theorem that replacing by does affect 
neither the consistency results nor the central limit theorem. Thus, whenever we refer to 
the variances of those random quantities, we use vi and i>2 without further notice. 

Before we start with the proofs of main results, we introduce more notation and state 
some simple lemmata. We consider the quantities 

= n^/H<Jim-i)/MWi^^ + Uif^), = n^^^<j^,^_,yMWj^% + i7,l^\), (6.1) 

which approximate ?m^'' and Y^^i, respectively, by using the associated increments of 
the underlying Brownian motion W. We further define 

frn _ l/iY{K) _ on pin _ 1/4t?(A:) _ r,ln ir, r,\ 

as the differences between the true quantities and their approximations. We further set 
/(x) := I a; I'' and g{x) := In the next lemma we study the stochastic order of the 
terms and i^^. 

Lemma 1. We have 

ma']+E[\c\']+E[\n'^'n^^n < c (6.3) 

for any q > 0, and 

E[\(3!^n + E[\f3':in + E[\n'/'Yi^^ ^ < C (6.4) 

for any Q < q<2{r + 1) + e with e as stated in Theorem 1. Both results hold uniformly 
in m. 

Proof. We begin with the proof of (6.3). In the case g > 1, we use Holder's inequality 
to obtain 
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■ E 



n/M -K+l 



nm/M-K 
i=n(j7i-l)/A/ 



< 



1 



n/M -K+l 

nm/M—K 

i=n(jri~l) /M 
1 

nm/M—K 
X ^ 

i=n(jri-l) /M 



11 



,1/4 



(i+;<-)/n 



i/n 



(i+_ff)/« 



ijn 



Thus the property follows from the boundedness of the functions a and cr, and a use of 
Burkholder's inequality. For g < 1, Jensen's inequality yields 

and we have (6.3) just as above. The corresponding assertions for and n^^^x'^^ can 
be shown analogously. 

Now let us prove (6.4). In the same way as before we have 

i?[|ni/^rWn < C(i?[|nV4f7(f)|^] +iJ[|„i/4x(f )|^]) 

for any g > 0. Moreover, it can be shown that n^^^U^'^ can be rewritten as a weighted 
sum of independent random variables, for which the convergence in distribution 

ni/4c7if)^AA(0,z.2C.2) 

holds. Using the continuity theorem and the moment assumption for each < q < 2{r + 
I) + e we obtain by uniform integrability of {ri^^'^u!^^ that £'[|ni/'*i7,tf'''|'] bounded. 
This proves (6.4) for n^/'^ym^' . The corresponding result for the quantities /3"j and /3^" 
can be shown analogously. □ 

The next lemma will be used later to prove that the error due to the approximation 
of ?Tr!^'' using is small compared to the rate of convergence. For a more general 
setting and a proof see Lemma 5.4 in Barndorff-Nielsen, Graversen, Jacod, Podolskij and 
Shephard [5]. 

Lemma 2. // 
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satisfies E[\Z';'J'] < C for all < q < 2{r + I) + e and if further 

M 

-J2E[\^^2~^,'Z'f]^o (6.5) 

m— 1 

holds, then we have 

1 ^ 



m — 1 

Proof of Theorem 1. Wc introduce the quantities 

M M 



MBV":=J2v?n and MBV'^ := J2 v'r^^, 



where t]^ and ry'" are defined by 

(r+l) /4 

C1C2 C1C2 ^ " 



Riemann integrabiUty yields 



so we are forced to prove 



and 



-L MB V"' -^MBViY,r,l), 



MBViY,r,l)n~ ^MBV" -^0 (6.6) 



-^{MBV - MBV'")^0 (6.7) 



in two steps. 

Considering the first step we recall the identity y/n = C1C2M and obtain therefore 



MBViY,r,l)n - j^MBV' = |J] (7^ - i;[7™|^(™-i)/M]), 

m— 1 

where 7m is given by 
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Using Lenglart's inequality and the ^(m_|_i)/M"™6asurablity of 7™ (for details see Lemma 
5.2 in Barndorff-Nielsen, Graversen, Jacod, Podolskij and Shephard [5]) we find that the 
stochastic convergence stated in (6.6) follows from 



J2E[\lm\^\^im-l)/M]^0- 

Holder's inequality and Lemma 1 yield i^[|7mP|-^(m-i)/M] ^ Cn~^, thus (6.6) follows. 

To prove the assertion in (6.7) recall that f{x) = |a;|'' and g{x) = \x\^ and observe that 
the continuity theorem implies 

£;[n(^+')/V(^(™-i)A/l^if' + uL^^)9i'^(m-i)/AM^^i + C/i+\)l^(™-i)A./] 

= MrMi(i'lO-^m-l)/A/ + '^2t^^)^''^'^^^ +Op(l), 

uniformly in m. Thus 

jriMBV- - MBV'n = V7 E i?[CI^(™^i)/M] + Op(l) 



M ' ' M 

m— 1 



with 



(r+l) /4 

C1C2 

To obtain the desired result it suffices to show 

1 

rn— 1 



which follows from 

M 

M 



M 

lirE^tlCI'l-O (6.8) 



m—1 

using the Cauchy-Schwarz inequality. In a first step we obtain for some constant C > 
ICP = i(/(C +/3™)5(C+i +/5™+i) - /(/?™)3(/3;"))' 
< C(.g'(C+i + /3;n+i)(/(C + PI.) - /(/?™))' 

+ /'(/^™)(5(c+i + - 5(/?:j.+i))' + fiPimP'^n+i) - 5(/3:;'))'), (6.9) 

where the quantities and are defined by (6.1) and (6.2), respectively. Since we have 
shown in (6.3) and (6.4) that the conditions on the boundedness of for an apphcation 
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of Lemma 2 are fulfilled, it suffices to prove 

M 
r/i— 1 

to obtain the assertion. 

For the first term in (6.10) we have 



(6.10) 



n/M -K + l 



n/M-K 

as in the proof of (6.3). Using (2.5) and 

{X{i+K)/n - ^tjn) - 0'(m-l)/M(W^(i+K)/n " W^i/n) 



(j+K)/n 



(i+A')/n 



Os ds + / (cTs - Cr(,„_i)/M) AWs 



i/n 



we obtain 



£'[|n^/"'((X(,;+A-)/„ - Xj/„) - 0-(„i_l)/A/(W^(i+A')/n - W^i/n))r 

r.(i+K)/ri 



i/n 



[cjs - cr(m-i)/M) ds 



(m-l)/A/ 



{(Ts - 0-(„,_i)/A/) di 



0(1). 



Consequently, 



= C^E I {as- a^Ms] /Aif d 



m— 1 



fo(l) 

0(1) 



C / i?[(a3-aLMsj/M) 1d.s + o(l) 



follows. Since a is bounded and cadlag, Lebesgue's theorem yields 

M 

M 



M 
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and as well for the second term in (6.10). We further have 

Since W^^\ is independent of at for any ^ < ^ we obtain 

m—1 rn— 1 

c 

- M E ^[l^™A-f - CTsI^ + - cr(m-l)/M| ]■ 
rn— 1 

The assertion therefore follows with the same arguments as above. That completes the 
proof of (6.7). □ 

Proof of Theorem 2. Theorem 2 can be proven by the same methods as Theorem 1. □ 

Proof of Theorem 4. Here we mainly use the same techniques as presented in 
Barndorff-Nielsen, Graversen, Jacod, Podolskij and Shephard [5] or Christensen and 
Podolskij [11]. We will state the proof of the key steps and refer to the articles quoted 
above for the details. 
We define the quantity 

M 

Lr^ir, I) = E (/(A" " E[f{(3l)g{P'„-J\T^,^_,yM]), (6.11) 

m—l 

where the terms and are given by (6.1), and f{x) — \x\^ , g{x) ~ \x\'- . In the next 
lemma we state the central limit theorem for Ln{r,l). 



Lemma 3. We have 



where L{r,l) is defined in Theorem 4- 
Proof. First, note that 

M+l 

L„(r,0= E C + Op(l), 

m=2 

where 0^ is given by 

C - "-i/4(/(/3;;_i)(.g(/3-_i) - M;('^i<-2)/a./ + ^^2^.')'^') 
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We have that 

'^[^ml-^(rn-l)/J\/] = 0, 

and 

A-f+l I r, o 2 2 rl 

CiC2 7o 



m=2 



V 



Next, let Z = W or Since 6*;^ is an even functional in W and B, and (M^, B) = -{W, B), 
we obtain the identity 

E[dmiZm/M — ■^(m-l)/M)|-^(m-l)/J\/] = 0. 

Finally, let N =^ (-^t)t6[o.i] be a bounded martingale on {il.,!F,{!Ft)t£io,i],P), which is 
orthogonal to W and B (i.e., with quadratic covariation [W^, 7V]t = [B,N]t = almost 
surely). By the arguments of Barndorff-Nielsen, Graversen, Jacod, Podolskij and Shep- 
hard [5] we have 

Now the stable convergence in Lemma 3 follows by Theorem IX 7.28 in Jacod and 
Shiryaev [20]. □ 

Now we are left to prove the convergence 

n^/*{MBV{Y, r, l)n - MBV{Y, r, 0) - i„(r, ^ 0. (6.12) 
Obviously, the convergence in (6.12) is equivalent to 

M 

P 



^^^0, (6.13) 

rn— 1 
M 



with 1?;;, 1?;'^ defined by 

= n-V4[/(^i/4y;W),,(^i/4y^W ) _ /(/3;jjg(/?-)], 

C - n^'' / {{-i-l + .2C.^)('-+')/^ - (-i^^™-i)/M + -20.'f^'~^'') An. 

J{ni-1)/M 

The convergence in (6.14) has been shown in Barndorff-Nielsen, Graversen, Jacod, Podol- 
skij and Shephard [5], and so we concentrate on proving (6.13). For the same reason as 
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in the proof of (6.6) (and from a similar argument as in the proof of (6.8)) this result 
follows from 

M 

^i?[C|.^(™-i)/A/]^0. (6.15) 

rn — 1 

Observe that 
Now we obtain 

i?[d-^(rn-l)/M] - E + C(2)|^(™-1)/m] + Op(l), (6.16) 

m— 1 m— 1 

with C»(l)> i^m(2) defined by 

C(l) = n-V4Vg(/3-)/(r.V4y W)e-, 

C(2)=n-i/4V/(/3;jJg(/3-)C, 

where are given by (6.2), and Wh denotes the first derivative of h. In fact, it is 

quite complicated to show (6.16) (especially when r or ^ G (0, 1]), but it can be proven 
exactly as in Barndorff-Nielsen, Graversen, Jacod, Podolskij and Shephard [5]. Note also 
that when r or Z G (0, 1] the terms V.g(/3^) and V/(/3^) are still well defined (almost 
surely), because CTs 7^ for all s. Assumption (V) implies the decomposition 

C = C(i) + C(2), 

where Cm(l)i Cto(2) are defined by 

^1/4 n/M~K ^(r+K)/n 

niM-K + l ^ (a„~a(™_i)/M)du 



i/n 

(i+K)/n / i^u 
ijn 



I a. As + / (cr^_ - cr[„_i)/A/) AWs 
J ijn J i/n 

^1/4 n/M-K . M+K)ln 



i=n{m-l)/M 



/ (K-V^;/„)dW„ , 

J i/n J 
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and a similar representation holds for Let us now prove that 

M 

^i?[i9:;(2)|^(„_i)/M]^0. (6.17) 

m— 1 

A straightforward application of Burkholder's inequality shows that 

M 

^ i?[V/(/3r„)5(/3;':)C(l)l-^(™-i)/A/] ^ 0. 

m— 1 

Next, note that since / is an even function V/ is odd. Consequently, V/(/3^).g(/3'")^," (2) 
is an odd functional of {W, V, B). Since {W, V, B) = -{W, V, B) we obatin 

M 
m — l 

which implies (6.17). Similarly we can show that 

M 

^£;[r^(i)|j-(™_i)/M]^o, 

m — 1 

which completes the proof of Theorem 4. □ 
Proof of Theorem 5. Theorem 5 can be proven by the same methods as Theorem 4. □ 
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